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We derive the steady-state electron distribution function for a semiconductor driven far from 
equilibrium by the inter-band photoexcitation assumed homogeneous over the nanoscale sample. 
Our analytical treatment is based on the generalization of a stochastic model known for a driven 
dissipative granular gas. The generalization is physically realizable in a semiconducting sample 
where electrons are injected into the conduction band by photoexcitation, and removed through 
the electron-hole recombination process at the bottom of the conduction band. Here the kinetics of 
the electron-electron and the electron-phonon (bath) scattering processes, as also the partitioning 
of the total energy in the inelastic collisions, are duly parametrized by certain rate constants. 
Our analytical results give the steady-state-energy distribution of the classical (non-degenerate) 
electron gas as function of the phonon (bath) temperature and the rates of injection (cw pump) 
and depletion (recombination). Interestingly, we obtain an accumulation of the electrons at the 
bottom of the conduction band in the form of a delta-function peak — a non-equilibrium classical 
analogue of condensation. Our model is specially appropriate to a disordered, indirect band-gap, 
polar semiconducting sample where energy is the only state label, and the electron-phonon coupling 
is strong while the recombination rate is slow. A possible mechanism for the dissipative inelastic 
collisions between the electrons is also suggested. 

PACS numbers: 05.70.Ln,72.20.jv,72.10.Di 

INTRODUCTION 

The kinetics of evolution of the electron distribution function for an electron-phonon system driven far from equi- 
librium by photo-excitation is of considerable current interest, experimentally [1] (the pump-probe experiments) and 
theoretically [2,3,4] (e.g., the two-temperature model. In the following, we develop a general stochastic approach 
to this kinetic problem as an alternative to the commonly used two-temperature model. Most such studies have 
been for metallic nanometric particles. More specifically, we derive the steady-state electron distribution function 
for the semiconducting sample driven far from equilibrium by the inter-band photoexcitation assumed homogeneous 
over the sample. Our analytical treatment is based on the generalization of a stochastic model known for a driven 
dissipative granular gas [5]. The generalization is physically realizable in a semiconducting sample where electrons 
are injected into the conduction band by photoexcitation and removed at the bottom of the conduction band through 
the electron-hole recombination process. Here the kinetics of the electron-electron and the electron-phonon (bath) 
scattering processes, as also the partitioning of the total energy in the inelastic collisions, are duly parametrized by 
certain rate constants. Our analytical results give the steady-state electron distribution function, and the mean energy 
of the classical non-equilibrium electron gas as function of the phonon (bath) temperature and the rates of injection 
(cw pump) and depletion (recombination). Interestingly, we obtain accumulation of the electrons at the bottom of 
the conduction band in the form of a delta-function peak — a non-equilibrium classical analogue of condensation! 
Our model is specially appropriate to a disordered, indirect band-gap, polar semiconducting sample where energy 
is the only state label, and the electron-phonon coupling is strong while the recombination rate is slow. A possible 
mechanism for the dissipative inelastic binary collisions between electrons is also suggested. 

THE MODEL 

Let ne{E)dE be the number of electrons lying in the energy range ztdE/2 centred about E in the conduction band 
of a semiconducting sample of volume Q. The electron-electron collisions, assumed inelastic in general, are described 
by the process; Ei + E[ — > Ef + E'^ = a{Ei + E'^) with a < 1, in which the tagged electron of energy Ei collides 
with another electron of energy E[ lying in the energy shell E[ ± jAE'^, and is scattered to the final state Ef. The 
scattering rate for this inelastic process is taken to be {1 — f)Tn{El)dE^. Similarly, the electron-phonon scattering rate 
is given by frnph{El)dEl, with nph{E[)dEl as the number of thermal phonons in the phonon-energy shell E'^ ± ^/^E[. 
Here, the fraction < / < 1 determines the relative strengths of the binary electron-electron and the electron-phonon 
collisions. Also, let the electrons be injected through photo-excitation into the conduction band at 
the rate gexS{E — Eg^), and then be removed(depleted) from the conduction band at energy Ed < E^x through 
recombination. Here the phonons are assumed to remain in thermal equilibrium at temperatures T. We expect this 
depletion-by-recombination to occur mostly at the bottom of the conduction band for which Ed ~ 0. This depletion 
rate can be modelled by a term —gdS{E — Ed)ne{Ed) . In our model sample we assume a uniform density of states for 
the electrons and the energy to be the only label for the single particle states. The photo-excitation is taken to be 
homogenous over the sample, which is reasonable for a nanoscale disordered semiconducting sample. For the above 
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dissipative model driven far from equilibrium, the kinetics for the non-equilibrium electron number density ne{E) is 
given by the rate equation 

= -neiE) J dE'[ne{E'){l - f) + nME')f]r 
+ dzp{z) J dE' j dE"S{E - za{E' + E"))ne{E')ne{E"){l - f)r 
+ dzp{z) j dE' j dE"6{E - z{E' + E"))ne{E')nph{E")fT 



In the above, we have assumed the total energy {E' + E") for a binary collision to be partitioned such that a 
fraction z, with probability density p{z), goes to the tagged electron of initial energy E', and 1 — 2; to the colliding 
particle (electron or phonon of initial energy E"). The inclusion of a in the electron-electron collision takes care of 
the possibility of inelastic electron-electron collisions. Note that we have suppressed the time argument {t) in the 
non-equilibrium electron- number density ne{E). Taking the energy Laplace transform 



heis) = / e-'^ne{E)dE, (2) 



of Eq. (1), we obtain, 



^n(s) = -rne(s)[(l-/)iVe + /iVj,ft] + (l-/)r^ p{z)dz nliazs) 

-\-fT I dz p{z) he{zs)hph{zs) 
Jo 

+geS)e-''''- - gae-'^'^neiEd). (3) 



In the following, we will consider for simplicity the steady-state condition under constant (cw) photoexcitation. 



gex{t) = gex- A pulsed excitation can, of course, be considered in general. Accordingly, we set M^he{s) = above 



and all quantities on the R.H.S. of Eq.(3) are then independent of time. 

In order to calculate the steady-state electron distribution for the system in terms of the bath (phonon) temperature 
and other rate parameters, we expand he{s) in powers of the Laplace-transform parameter s as 

ne{s) = N,-s{Ee)+s\El)..., (4) 

and equate the co-efficients of like powers of s. Thus, from the zeroth power of s, we obtain at once 

ne{Ed) = {gex/gd}- (5) 

Similarly, from the first power of s, we get. 



, s. ^ ifmjeph) .9e.(ge. -i?rf)/r 

Pe-phil-aKl-f) + f/2 7V2^Pe-p4Pe-ph(l-a)(l-/) + //2]- ^' 

In the above, we have taken a uniform limit for the energy partition: p(z) = 1. 

Here, we have defined (ce) = {Ee)/Ne = mean electron energy; {Cph} = {Eph)/Nph = mean phonon energy 
(= UbTb); and Pe-ph = Nf,/Nph = electron-to- phonon number ratio. It is to be noted that in the limit a = 1 {i.e., 
for elastic electron-electron collisions as is usually expected for an electronic system unlike the case of the granular 
gas), and g^x = {i.e., no photo-excitation), we recover (eg) = {epu), i.e., the electrons and the phonons are at the 
same temperature, as is physically expected under equilibrium conditions. In general, however, the mean electron 
energy in the steady state is not the same as the mean phonon energy, and the former depends on the excitation rate 
(the drive gex)- In the following, we will consider other choices for p{z). 

Case 1 {Extreme Partition Limit): p{z) = \{5{z) -\-5{z — 1)) 
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Analytic form for ne {s) in the steady state can be obtained for the particular case of clastic scattering, {a = 1), and 
for the partitioning function p{z) = ^S{z) + ■^S^z — 1). It can be easily verified that for this insertion of the Eqs. 
(5) and (6) remain unchanged. So, with this choice, Eq. (3) becomes 



he{s)[{l - f)N, + fN,t,] = (!-/) + \5{z - l))dzhl{zs) + / 

1 

*-®<*]/r. 



+ - l))dzne{zs)nph{zs) + gex[e 



— e 



This can be readily solved to give. 



with 



2(1-/) 



C = {1- f)Nl + fN.Nph + 2fc(e- 



2(1-/) 



(7) 



(8) 



)/r. 



Wc have to choose the —ve sign in the above solution so as to satisfy ne(0) = N^, nph{0) = Nph- Here, we have 
assumed that the bath phonons obey the equilibrium Boltzmann distribution {i.e., nph{s) = Nph/{1 + {Eph)s)). It is 
to be noted here that he{s) is not zero in the limit s — > oo. This is due to the fact that there is a delta-function peak 
of strength 



V = Nphi 



-ph + 



/ 



1-/ 



-ph + 



f 



f 



at energy E = 0. The strength of the delta function reduces to Nph 



(9) 



when there are equal number of electrons 



and phonons in the sample. It is to be noted that we have here three macroscopic unknown quantities, Ne,< Eg >, 
and the strength rj of the delta-peak at i? = 0, while there are only two equations namely, Eq. (6) and Eq. (9). It is, 
therefore, necessary to solve Eq.(8) for the full distribution fie(s) in order to determine these macroscopic quantities. 
We have, however, a simpler case when the depletion by recombination is at the bottom of the band, i.e., Ed = 0. 
In this case, the last term in the rate equation (1) becomes —gd5{E)ne{E). Thus, for a delta-function peak in nf,{E) 
of strength rj a.t E = 0, this term becomes —rjgd5^{E). Recalling that S'^{E) = £^5{E)/2tt with ^ an infinitely large 
conjugate interval, we can redefine a renormalized depletion rate constant gd^l^i^ = Gd-, giving the depletion term as 
—r]GdS{E). This enables us to determine the condensate strength r] directly in terms of the photo-excitation rate gex- 



Case II (Equipartition Limit): p{z) = 5{z — ^) 

Finally, we solve Eq.(3) for the steady-state condition for elastic collisions (a = 1) with equipartition of energy 
(z = 1/2) between the colliding particles. In this case also, it can be easily verified that with this insertion of .2 = 1/2, 
the Eqs. (5) and (6) remain unchanged. For this special case, Eq.(3) reduces to 



Now, in the limit of small s, the he{s) can be written as, 



ne(s)[(l - f)N, + fNph] = (1 - f)nl{s/2) + fn,{s/2)nph{s/2) + ^(e'^^- - e-^^-). (10) 



poo nOC 

ne{s)= / e-'^ne{E)dE / {1 - {s/2)E...)e-^'/^^^ne{E)dE 
Jo Jo 

c^n,{s/2)-{s/2){E,)-{s/2f{El). (11) 
Using Eqs. (10) and (11), we obtain the electron distribution function in the s-domain as. 



"«^^/2)- 2(1-/) + 2(1-/) ^^'^ 

with 

Ci = s{E,)[{l - f)N, + fNph]/2 + fle.(e-^^'^ - e-^^^)/r 
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To invert the above expression for the electron-number density distribution from the s-domain to the i?-domain, 
we have carried out the small-s analysis. Equation (12) can then be written as 



ne{s/2) 



2Ne 

{Eph) 



1 



1 



(Eph) \s + 2/{Eph) 



[(1 - f)N, + fNph]{Eph)N, - 2fci 



il-f){Eph)N, 



(13) 



with 



We finally obtain 



ne{E) ~ 
X ls{E) 



N, 



ki = {Ee)[(l - f)N, + fNph]/2 + gex{Eex - Ed)/r. 

(Ee) 



f 

1-/ 



Nph - 



(Eph) 



{Eph) 
AN, 



f 



2gexiEex — Ed) 



{l-f)p,^phj [1- f)T{Eph)N, 



(Eph) 



(14) 



Clearly, it has two parts the Boltzmann distribution at temperature T , and a delta-function condensate at energy 

E = Q. 



DISCUSSION 

We have treated here the problem of energy-distribution of photo excited electrons in a semiconducting sample as 
a generalization of the model for dissipative granular gas driven far from equilibrium. An interesting feature of the 
above non-equilibrium distribution of the cw photo-excited electrons is the delta-function peak (classical condensate) 
appearing at the bottom of the conduction band. It should be possible to probe this steady-state feature through a 
pulsed pump-probe experiment combined with the cw excitation. Of course, for the full electron-distribution function, 
the Laplace transformed ne(s) has to be inverted numerically. It is to be noted that we have assumed the sample to be 
photo-excited homogeneously {i.e., only spectral but no spatial diffusion). This is clearly relevant to nanometric-scale 
samples. As for the dissipative feature (a < 1) of the binary electron-electron collisions, it can possibly derive from 
the dissipative polarization of the dielectric medium that mediates the electron-electron Coulombic interaction [6]. 
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